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Abstract 

Employing probabilistic techniques we compute best possible upper and lower 
bounds on the price of an option on one or two assets with continuous piecewise 
linear payoff function based on prices of simple call options of possibly distinct 
maturities and the no-arbitrage condition, but without any assumption on the price 
dynamics of underlying assets. We show that the problem reduces to solving linear 
optimization problems that we explicitly characterize. We report numerical results 
that illustrate the effectiveness of the algorithms we develop. 
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1 Introduction 



One of the central questions in financial economics is to find the price of a derivative 
security given information on the underlying asset. Pricing of a derivative is usually 
realized by solving the Black-Scholes equation [3], which is based on the assumptions 
of a Geometric Brownian motion for the price of the underlying asset and no-arbitrage 
in the market. One of the parameters in the equation, the volatility of returns on the 
underlying asset, is assumed to be defined and constant, while in reality it is varying in 
time and not known in advance. Practitioners typically use the implied volatility to value 
a derivative security by the Black-Scholes formula, that is they extract the volatility from 
prices of other options in the market, based on the assumption that the Black-Scholes 
equation is valid and other options are correctly priced. 

A natural question that arises is to determine the range of values for the price of an 
option based only on prices of other options and the no-arbitrage assumption, but with 
no assumption on the price dynamics of the underlying asset. 

In their seminal work, Cox and Ross [7] and Harrison and Kreps p~2] show that the 
condition of no-arbitrage is equivalent to the existence a probability measure Q, equiva- 
lent to the original measure P, with respect to which all discounted securities processes 
are martingales. Rubinstein [TS] and Longstaff [16] introduce the idea of deducing the 
martingale measure from observed European call prices by solving a quadratic optimiza- 
tion problem. 

Lo [TS] derives best possible closed form bounds on the price of a European call option 
given the mean and variance of the underlying stock price under risk neutral measure. 
Grundy [UJ extends Lo's work for the case when the first and the kth moments of the 
stock price are known. Bertsimas and Popescu [3] derive best possible bounds of the 
price of a European call option, as well as on moments of the prices of the asset, given 
prices of other similar options, using a convex optimization approach. D'Aspremont and 
El Ghaoui p] address the problem of computing upper and lower bounds on the price 
of a European basket call option, given prices on other similar baskets. They introduce 
a linear programming relaxation of the problem and show that this relaxation is best 
possible in some special cases. All the above problems are solved for the case of a single 
maturity and single underlying stock. 

In this paper, we consider the problem of determining upper and lower bounds on 
the price of a European option given prices of other options of different maturities. 
In particular, the contributions of the article are the following: 

(a) We establish necessary and sufficient conditions that European call options of dif- 
ferent maturities should jointly satisfy, so that there exists no arbitrage. From this 
structural property and using geometric arguments, one can derive best bounds on 
the price of a European call option with a given strike price and maturity. As we 
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present more general methods in later sections we do not present this method here 
(see Bertsimas and Bushueva [2]). 

(b) Given a set of European call options of different strike prices and different maturi- 
ties we determine best possible upper and lower bounds on the price of an option 
with a continuous piecewise linear payoff function. Based on direct probabilistic 
methods, we solve the problem as a linear optimization problem. 

(c) Given European call options on two individual assets we determine best possible 
upper and lower bounds on the price of an option with a payoff being a continuous 
piecewise linear function of prices of the two assets at the time of option's maturity. 
Again, options are allowed to have different maturities. We present two algorithms 
based on linear optimization. The first is an exact formulation and involves a linear 
optimization problem of potentially exponential size; the second is asymptotically 
exact and involves a linear optimization problem of polynomial size. 

(d) We determine upper and lower bounds on the price of a European basket call 
option, given prices on other similar baskets solving exactly the problem addressed 
in D'Aspremont and El Ghaoui 1 approximately. 

Throughout the paper, we refer to these problems as Problems (a), (b), (c), (d). The 
paper is structured as follows: in Sections 2, 3, 4, 5 we present our solutions for Problems 
(a), (b), (c) and (d) respectively. 



2 Characterization theorem for the case of multiple 
maturities 

In this section, we determine best possible bounds on the price of a European call option 
given prices of options on the same underlying security, but with potentially different 
maturities. 



The price of a European call option with a strike price k and maturity t is given by 

/ 



C = 



r(s)ds 



(S t - ky 



Eq 



r(s)ds 



V 



(1) 



where r(s) is the instantaneous riskless rate of return, {St}t>o is the price process of 
the underlying security defined on some probability space (f2, B, P), and Q is a measure 
equivalent to P, such that {A t } t > := I S t e~ So r ( s ) ds I j s a martingale process under 
Q. Two measures are equivalent if they have the same null sets. 
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Thus, the problem we address in this section can be reformulated as follows: given 
the set 

S := {(fcy.Cy,^) 6l 3 + |*i <t 2 <... < t n ,i€{l,2,...,n}, j € {1, 2, . . . , U(i)}} 

find the sufficient and necessary condition that values of the set should jointly satisfy so 
that there exists a nonnegative martingale {X tl , X t2 , . . . , X t > , . . . , Xt n }, such that 



E 



— cv 



(2) 



for all points of the set. 
Remarks: 



(a) We can set without loss of generality r(s) = 0, since we can define kij — h^e ^ r i s ) ds : 



and rewrite Eq. ^ as E 



(x u -k 3 y 



Ci. 



(b) Since we can assume r(s) = 0, only the order of t%, t2, ■ ■ ■ ,t n is important, but not 
the specific values that the sequence takes. 

(c) The price of an option with a strike price is given by E[X t ] = X , which is the 
current stock price and, thus, is always known. 

In light of the above remarks, a simpler formulation of the problem is as follows. 

Problem (a) 

Given an ordered set 

S := {(k tj , C tj ) e K 2 + | t e {1, 2, . . . , n}, j S {1, 2, . . . , U(t)}} (3) 

and Xq e M + find the sufficient and necessary condition that values of the set should 
jointly satisfy so that there exists a nonnegative martingale {X\, . . . , Xf , . . . , X n } such 
that 



E 



(X t - k tj ) + = C ts for all (ktj, C tj ) e 5, 



(4) 



and 



E[X t ]=X Q for alU= 1,...,; 



Let C designate the class of probability laws with support in R+ and finite expecta- 
tions. We associate with each element II of C the transform : R+ — > R+, defined 

as 

/•OO 

*n(i) = E n [(x - t)+] = (x- t)+dU(x). (5) 

Jo 

We next show that ^u(t) uniquely determines II. 
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Theorem 1. (a) The right derivative * n '(*+) of \l>n(*) exists for alH > and 

*n'(t+) = -H( (t,+oo] ). 

(b) Let III, n 2 S C and for each f e M+, ^(i) = *n 2 (*)- Then, IIj = n 2 . 
Proof, (a) Let t > and ft > 0. We have 

^ = J ^ dU{x)=J Jt,h( x ) dU(x), 

where the function ft.h{ x ) is defined as follows (see also Figure[lJ: 

0, x <t, 

(x-(t + h))+-(x-t)+ I x-t 
ftA x )~- " l " —={ —> t<x<t + h, 

-1, x>t + h. 



ft,h( x ) 


t t + h 




\ 


X 


-1- 









Figure 1: The function ft,h{ x )- 



Since the sequence of functions ft,h( x ) is monotone (in h), and ft,h( x ) converges as 
h — > + , then by by the monotone convergence theorem [51 p. 100], J + °° ft,h(x)dH(x) —> 
-n((i,+oo)) as h -> 0+. Thus V n (t+) exists and equals — +oo)). 
(b) Immediately follows from (a). □ 

The following proposition summarizes the properties of the 'J'— transform. 

Proposition 1. Let II £ C. Then ^n(t) satisfies: 

(a) ^nit) is a nonincreasing function of t S M. + ; 

(b) ^n(i) is convex; 

(c) lim^+oo ^u(t) = 0. In particular, if IT has a bounded support, then ^u(t) = 
for t large; 

(d) * n (t) > * n (0)-t. 
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Proof, (a) Notice, that ^u(t) is a continuous function: 



|*n(* + fc)-*n(*)l = 



(x - t)dTI(x) 



<\h\. 



Since \&n(£) is continuous and \f r n'(*+) = — H( (t, +00] ) < 0, then \I/n(i) is nonincreas- 
ing. 

(b) /^(i) := (x — t) + is convex in i for each x. Thus, *£n(t) is convex as a convex 
combination of convex functions. 

(c) 

*n(i) = / {x -t)+dU(x) = / (x - t) + dH(x) < / xdU(x). 
Jo Jt Jt 

Since Jq + °° xdH(x) < 00, then J t + °° xcZII(x) — > as i — > 00. Thus, ^n(t) — > as t — > +00. 
(d) * n (t) = .CO - *)+dn(x) > f °°(x - f)dTI(x) = # n (0) - 1. □ 

We next show that the properties in Proposition [1] are also sufficient for a function 
to be a 'J— transform of some distribution II. 

Theorem 2. Let g : R+ — > K + be a function such that (a) 3(2) is nonincreasing, (b) 
(7(2) is convex, (c) lim a; _ ( + 00 g(x) = 0, and (d) g(x) > g(0) — x. Then, there exists II € C 
such that g(x) = ^u(x). 

Proof. From the convexity of g it follows that the right derivative of g, g'(x+), is a right 
continuous nondecreasing function. Since g(x) > g(0) — x and g (x) is nonincreasing, then 
— 1 < .9'( a; +) < 0- Also we have <7'(x+) — * as x — > +00. To justify the last statement, 
notice that </(x+) is a monotone bounded function, and consequently, it must converge 
to some finite value a as x — > 00. However, o ^ would contradict the requirement in 

(c) . 

Consider the function F(x) := l + g'(x+) defined on R + . From the above statements, 
it follows that (i) < F(x) < 1, (ii) F(x) is a right-continuous and nondecreasing, and 
(Hi) F(x) — > 1 as x — > 00. Thus, -F(x) satisfies all the properties of a distribution 
function. Denote by II the law with the distribution function F(x). Then, 

r + co r + co 

*n(0)=En[X}= (1-F(x))dx= (-g'(x+)) dx = g(0) - lim g(x) = g(0). 

Jo Jo 

The second to last equality (identified with "*") holds since g(x) is a convex function 
and thus its right and left derivatives might differ only on a countable set of points. 

Since $n'(^+) = — (1 — F( x )) = 9'( x +): ^n(^) and g(x) are both continuous and 
* n (0) = 3(0), it follows that ^u(x) = g(x) for any x S M+. □ 

Since there is a one to one correspondence between laws of random variables and their 
'J— transforms, it follows that conditions that a sequence of laws of random variables 
should satisfy can always be reformulated in terms of conditions on the sequence of 
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transforms of these random variables. Theorem [3] establishes necessary and sufficient 
conditions that the sequence of 'J— transforms should satisfy so that the corresponding 
random variables form a martingale. However prior to proving that theorem, we will 
need the following auxiliary results. 

Definition 1. Let III, II2 <E C. Define the relation IIi7?.Il2 if and only if there exist 
random variables X\ and X2, defined on the same probability space, such that {X\, X2} 
is a martingale and X\ has a marginal law III, while X2 has a marginal law II2. 

Proposition 2. The following properties hold: 

(a) If IIiftn 2 and TI2IUI1 then III = I± 2 ; 

(b) ni^Hi; 

(c) if ni^n 2 and n 2 rai 3 , then n"i7en 3 . 

Proof, (a) If {Xi,Bi}, {X2,B 2 } is a martingale with marginal laws IIi and II2 corre- 
spondingly, and / is a convex function, defined on R, then f(Xi) = f(E [X2\Bi]) < 
E [f(X2)\Bi] by the conditional Jensen's inequality, and thus E[f(Xi)} < E[f(X2)]. 
(X — 1) + is a convex function of X, so ^rii (t) < \Pn 2 (*) f° r each t £ R + . From Il27?.IIi by 
the same argument we obtain 1 $n 2 (t) < ^rii(i) for each t £ M + . Thus ^riiC*) = ^n 2 (*) 
and by Theorem [TJ LTi = II2 . 

(b) {Xi,Xi} is always a martingale. 

(c) Let {Xi, X2} be a martingale with the law C(Xi, X2) defined onlxl with support 
in R + x R + and marginal laws IL and II2 respectively. Let {X2,X 3 } be a martingale 
with the law C(X2, X3) defined on M x M with support in R + x R + and marginal laws 
1T2 and II3. Then on R there exist conditional distributions C(Xi\X 2 ) for £(X 1; X2) 
and £(^31X2) for £(X2, X3). Since the marginal laws of X2 and X2 are equal, then 
by Vorob'ev-Berkes-Philipp theorem [10, p. 7], we can define a law £(Xi, X2, X3) on 
R x R x R, such that X\ and X3 are conditionally independent given X2, that is 

L{{X X ,X Z )\X 2 ) = C{Xx\X 2 ) x C(X 3 \X 2 ). 

Denote by Bj, i = 1, 2, 3, the smallest a-algebra for which all -Xj., k < i, are measurable. 
Then {Xj,Bj}? =1 is a martingale with marginal laws IIi,Il2 and II3 correspondingly. 
Note that the martingale thus constructed is also a Markov process. □ 

Theorem 3. Let Hi, LI 2 , . . . , II„ £ C and E n , [X] = E Ul [X] for each i = 1, . . . , n. Then, 
there exists a martingale {Xi, X2, . . . , X n }, such that Xi has marginal law IIi for each i 
if and only if the sequence of functions ^xi t > *n 2 , ■ ■ ■ , is nondecreasing. 

Remark: The theorem was obtained in Kertz and Rosier [14 . The key step in their 
proof is a result of Strassen (T^]. Our approach is completely different. 



7 



Proof. The only if part is a direct consequence of the fact that (X — t) + is a convex 
function of X and Jensen's inequality for convex functions of martingales [9j p. 277]. The 
proof of this statement is also in the proof of part (a) of Proposition O 

In the other direction, we have to show that it is possible to define a joint distribution 
of X\, X2,- ■ ., X n with the given marginal laws (since they are uniquely determined by 
*rii) ^nsj ■ • ■ 7 ^n„) ; so that this sequence is a martingale. It is enough to show how to 
define a joint distribution of X\ and X2, since afterwards we can proceed recursively to 
form a martingale which is a Markov process as in the proof of Proposition [21 

We first define a sequence Xi, Yi, Y2, . . ■ , Y n , . . ., so that this sequence is a martingale 
and \fV n — > ^n 2 ■ We will then show that there exists Yoo, such that {Y n }i<n<oo is a 
martingale and Yx, has a distribution II2 . Then X2 can be defined as X2 = Yx ■ 

Let Y"o = X\ and let us form a martingale Yq, Yi, Y2, . . . , Y n such that it is also a 
Markov process, and thus, we only need to define the joint distribution of Yi and l^+i 
for each i. Now, the joint distribution of Yi and Yi+i is uniquely determined as soon 
as we know for each y £ K+ the conditional distribution of Y + i given Yi — y. Thus, 
we are going to define the conditional distribution for Y\, then for Y2 and so on. We 




a b V T t 

Figure 2: Construction of 4^ . 



know that — < ^nv Now Y\ will have its ^-transform, equal to 
everywhere except on (a, b) where it will be a straight line (Figure^) We choose (a, 6) so 
that the line which goes through the points (a, (a)) , (b, (b)) is below the function 
^1X2 (*) (the specific choice of (a, 6) is discussed later.) What distribution will have its 
^-transform equal to the transform of IT everywhere except (a, b) and being a straight 
line on (a, 6)? From TheoremQ] (a), we see that should have no weight on (a, b) but 
the weight of [a, b] should be the same as for the distribution III. Thus, we take all the 
weight of (a, b) and redistribute it between {a} and {&}. We have to do it in such a way 
that {Xi, Yi} is a martingale. 

The law C(Xi) is uniquely determined by ^ x t and the law C{Y\) is uniquely de- 
termined by ^yi- We define the law P := C{X\,Y\) with support in R + x M + in the 
following way. As the marginal law of X\ take • Then for each X\ £ R + define the 
conditional distribution of Y\ given X\. First, if X\ £ {R + \(a, &)}, take Y\ = X%. Now, 
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if Xi £ (a, b) we should have P(Y\ € {a}U{6}|Xi £ (a, b)) = 1 for (t) to be a straight 
line on (a, b) and equal to 'J'xi(i) everywhere else. Let us denote 7rx 1 (^i,5 / i) := X\. 

Thus, for any c £ {a,b) we want the joint distribution of X\,Y\ to be such that the 
following holds 

aP{Y 1 = a|Xi = c) + bP{Y 1 = b\X x = c) = c, 
P(Yi = a|X x = c) + P(Yx = = c) = 1. 

Consequently, 

P(Y 1 = a\X, = c) = h -^- , P(Y 1 = b\Xi = c) = C --^. 

b — a b — a 

Thus, the conditional distribution of Y\ given X\ is defined and {X\, Yi} is a martingale. 

In the same manner we can define the conditional distribution of li+i given Yj, for 

i = 1,2,.... 

The next step is to show that \& x 2 (*) can be approached by ^y„ (*) with any given 
precision, that is that we can choose Y l5 Y 2 , . . . , Y„ so that ^y„(i) ~ * ^x 2 (*)- For each 
i, 1 i>Y i+1 {t) is different from \Py;(t) only on a certain interval (ai,bi), where ^Y i+1 (t) is 
a straight line. Thus, the choice of the sequence {Y n } is equivalent to the choice of the 
sequence {(a„, &„)}■ 




ai t\ bi 

Figure 3: The choice of (ai, 6i 



Let 

[/:={feR + :$x 2 (t)>%(t)}. 

Let be a countable dense set in U. At each ij, i = 1, 2, . . . , draw a tangent line j/i(t) 
to *x 2 (0- ( See Figure^) Let us first assume that *x 2 (*i) 7^ °- Then #x- 2 '(**+) < 
and the tangent line intersects Vtyj^ft) at two different points, <ij € M+ and &i € R+, 
since (i) -> as t -> and ^y (0) = 1, while is below <fx 2 (*) and thus j/j(0) < 1. 
Thus, the choice of {£„} is equivalent to the choice of {(a n ,b n )}. We define = 
max{j/i(t),*y 4 (t)}. In case *x 2 (*i) = 0, \&x 2 + (*») = and we take := 

Thus, the sequence {^y„} is an increasing sequence of continuous functions converg- 
ing to a continuous function ^x 2 on a countable dense set. Consequently, \ty n — > f x 2 
pointwise on R+. 
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Since ^x x (t) — > and ^x 2 ^ as t ^ +00, then for any e there exist T £ M. + 
such that for any t > T, \^x 2 (t) ~ ^Xi(t)\ < e - Then for any n and any t > T, 
\^x 2 (t) — ^Y„(t)\ < e - Taking into account that {^y,,} converges uniformly to 1 $>x 2 011 
[0, T], we conclude that \Ey n converges uniformly to ^ x 2 on R + . 

Now, since — > ^ x 2 uniformly on R + and the right derivatives of $>Y n (t), n = 
1,2, ... , and $ x 2 are monotone functions, then also VPy n / (i+) — > , J r x 2 (*+)• Thus, the 
distribution functions of Y n , n = 1, 2, . . ., converge pointwise to the distribution function 
of X 2 . 

We next show that {Y n } is uniformly integrable, i.e., for any e > there exists to > 0, 
such that for all t > t , 

supE [Y n l {Yn>t} ] < e. 

n 

For each t and n we have ^>Y n (t) < 1 &x 2 (t)- Therefore, 

r+<x r+00 

/ (Y n - t)+dU Yn < (X 2 - t)+dU X2 , 
Jo Jo 

/+00 r+oc 
Y n dn Yn < J x 2 du X2 +t\u Yn ( [t, +00] )-n X2 ( [t, +00] ) 1 . 

Since X 2 dUx 2 — ^ as i — > +00, then we can take to such that / t + °° X 2 dllx 2 < e/3 
for all t > to- Now since H Yn ([to, +00]) — > llx 2 ([ioj +00]), there exists no such that for 
all n > no, |IIy n ({to, +00]) — TLx 2 ([to, < gf - - Thus, for any n > no and t > to we 

have 

£ [^V„> t} ] < £ [^l{r„>t o} ] < I + = I ■ 

Therefore, 

r i 2e 

sup E [Y n l {Yn>t }\ < — • 

n>no «J 

For Yi, I2, . . . , Yn, we can always choose t\ big enough so that for all t > ti , E [Yils Yi >t}\ 
< |, z = 1,2,... ,no, since all Yi have finite expectations. Consequently, if we denote 
t* := max{to,£i}, then 

supE [Y n l {Yn>t} ] < e, 

71 

for all t>t*. 

Thus, Yi,Y 2 , . . . ,Y n is uniformly integrable, and, consequently [H p. 283], right clos- 
able, which implies that there exists a random variable Yoo such that _E[Yoo|6„] = Y n for 
all n, where B n designates the smallest a-algebra for which Y\, . . . ,Y n are all measurable. 
Then by Doob's theorem |J} p. 285] Y n converges to Y x a.s., and consequently the distri- 
bution functions of Y\ , Y 2 , . . . , Y n also converge to the distribution function of Yoo . But 
then Yoo and X 2 have the same distribution function, and thus the same ^-transform. 
So we take X 2 = Y^ and X\, X 2 is a martingale. 

□ 
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Now we are ready to formulate the conditions that the set S defined in ([3]) should 
satisfy, so that there exists a martingale {X\, X2, • • ■ , X n }, such that Eq. (j4|) holds. 

Proposition 3. The set 5, as defined in ((3]), satisfies the condition of martingale exis- 
tence if and only if there exists a nondecreasing sequence of convex functions, {fft(a;)}tg{i.... ! n} : 
K + — * R + , such that for each t £ {1, 2, . . . , n} 

(a) 9t (0) = gi(0); 

(b) g t (x) > g t (0) - x; 

(c) gt(x) — > as x — > 00; 

(d) g t {k tJ ) = C tj for each j £ {1, 2, ... , £/(*)}. 

Proof. " " By Theorem [2] for each g t there exists II ( S C such that ^n t = 9t- Then, 
since the sequence of functions ^rii , ^ru , • • • , ^n„ is nondecreasing and \&n t (0) = ^rii (0) 
for each t, from Theorem [3] it follows that the martingale satisfying all the points of the 
set S exists. 

" " Suppose that there exists a martingale {Xi,X2, ■ ■ ■ ,X n }, such that holds. 
Then we can take gt(x) := ^x t (x) for t = 1, 2, . . . , n. By Theorem [3] the sequence of 
functions {gt}t=i,....n is nondecreasing and by Proposition [1] conditions (a)-(c) hold. By 
definition, g t (k t j) = * x t {hj) = C t j for all (k tj , C t j) £ S. □ 

Designate by 

Sf.= { (klj , Qj ) I (fcy , Cy ) G 5, / = *} (6) 

the set of options of maturity t and by 

S> t := U;> t «S ; (7) 

the set of options of maturities t and higher. Let also 

5g :=5> t U{(0,+co)}U{(+co,0)}. (8) 

Finally, let <96>>t : R+ — > R+ represent the lower boundary of the convex hull of 
U {(0,X )}, that is 

dS> t {x) = mm{y\(x,y) £ convex hull of U (0,X )}} . (9) 

Theorem 4. The set S satisfies the no-arbitrage condition if and only if 

(a) C tj > X Q - k tj for each (ktj, C t j) £ <S; 

(b) for each t, none of the points of St is in the interior of the convex hull of S> t . 

Proof. " => " Define functions <?t(fc) : R+ — > R+ for each t = l,...,n, such that 
gt(k) = dS>t(k) Vk > 0. Then the sequence of functions {gt} satisfies all the properties 
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of Proposition [3] except of property (c). Each function g t is piecewise linear with slopes 
of all, but the last, line segments being negative. The last line segment has a slope zero. 
Let Xt := min{a; | x € K+, g't(x) — 0}. Denote by s the set of all slopes of functions gt, 
t = 1, 2, . . . , n, and let s* := max {s\s £ s, s < 0}. Define the function 



9t(x), 

gt(x) := { g t (x t )+s*(x-x t ), 

0, 



x < X t , 

{x > x t | g t (x t ) + s*(x - x t ) > 0} , 
otherwise. 



Then the sequence of function gt{x) satisfies all the properties of Proposition[3]and, thus, 
the martingale exists. 

" <== " If Xi, X%, . . . ,X n is a martingale, satisfying all the points of S, then for any 
(k, C) which belongs to the convex hull of S>i, it must be that C > , i'x 1 (k) > Xq — 
k by Proposition [1] and Theorem [3] Assume that there exists a point (kfj*, Ct*j*), 
t* £ {1, . . . , n}, j* £ 1, . . . , U (t*), which is in the interior of the convex hull of S> t , 
(Figure 12) Take points (k l ,C l ) £ S^ t , {k r ,C r ) £ S? t such that dS> t *(k l ) = C l and 
dS>t*(k r ) — C T and k l is the closest to kfj- from the left, while k r is the closest to kfj* 
from the right. Then {kfj*,Cfj») lies strictly above the interval connecting (k l ,C l ) 
and (k r ,C r ). Since V Xt » < *x t , +1 < ■■• < by Theorem^ then * Xt , (k l ) < C l 
and ^x t * (k r ) < C r , thus (kfj*,Cfj») is also strictly above the interval connecting 
(fc^^Xj, (k )) and (k r ,^x t * {k r )), which contradicts to the condition that *$>x* is a 
convex function. □ 



C 




Figure 4: The point {kfj*, C t *j*), t* £ {1, . . . , n}, j* £ 1, . . . . U (t*) is in the interior of 



the convex hull of S> t , . 



From the last theorem, we can find bounds on the price of a call option, given the set of 
priced options based on a geometric approach, see Bertsimas and Bushueva [2]. However, 
we will suggest another approach, which leads to the linear optimization problem, and 
that is extendable to multiple dimensions. 
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Restriction to Markov martingales 

The idea of our approach is to reduce the class of martingale distributions over which we 
maximize (or minimize) as much as possible, while making sure that the optimal solution 
does not change. In this section, we show that it is enough to consider martingales with 
the Markov property only. 

Proposition 4. Let X\,X%, . . . , X n be a martingale with values in M. d and let \ii be the 
marginal distribution of . Then, there exists a martingale with the Markov property 
X\, X%, • • • , X n , such that the marginal distribution of Xi is \ii for all i = 1, 2, . . . ,n. 

Proof. For specifying the distribution of a discrete time Markov process, it is enough to 
specify the transition probabilities as well as the initial distribution. 

Let v = C(Xi), fXi t 2 = C(X2\Xi), . . . , /i n _i,„ = C(X n \X n -x). Then by the existence 
theorem for Markov processes due to Kolmogorov [13l p. 120] there exists a Markov 
process Xi, ...,X n such that C{X X ) = v, C{X 2 \X X ) = /i 1)2 , C(X n -x\X n ) = /z„_i, n . 

It remains to show that the so-defined process X\,X 2 , ■ ■ ■ ,X n is also a martingale. 
Let a k be the smallest a-algebra for which X\, X 2 , ■ ■ ■ , X k are all measurable. We need 
to prove that E[Xk+i \<Jk] — Xk- From the Markov property of Xi, . . . ,Xk, it follows 
that E[X k+1 \a k ] = E[X k+1 \X k ]. By definition E[X k+1 \X k ) - E[X k+1 \X k }. 

On the other hand, since X%, X 2 , . . . , X n is a martingale, we have 
E [E[X k+1 \X k ,X k ^, . . .,Xx}\X k ] - E[X k \X k ] = X k . Also 

E[E[X k+1 \X k ,X k _ u ...,Xi]\X k ] = E[X k+1 \X k }. Thus E[X k+1 \X k ] = X k . Conse- 
quently, E[X k+1 \X k ] = X k and E[X k+1 \a k ] = X k . □ 

Definition 2. A stochastic process which is a martingale and a Markov process will be 
called a Markov martingale. 

Corollary 1. If we restrict our consideration to Markov martingales only, the optimal 
solution to our problems does not change. 

3 The One Dimensional Case 

In this section, we address Problem (b), i.e., we find best possible upper and lower 
bounds on the price of a European style option with a continuous piecewise linear payoff 
function, given the set of European call options of different maturities. 

The price of a European call option with a strike price k and maturity t is given by 
Eq. Ilj. The problem we address in this section can be formulated as follows: given a 
continuous piecewise linear function g : R+ — * R+ and a set 

S := {(kij, Cij, t{) | {kij, Cij, ti} £ R + , ti < t 2 . . . < t n , 
i £ {1,2,. ..,n}, j£ {1,2,..., [/(i)}} 
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find the domain of values of 



e Jo 



r(s)ds 



E 



under the condition, that {X tl , X t2 , . . . , X t » , . . . , X tn } is a nonnegative martingale such 
that 



E 



X ti -kije Jo 



u . ■ - 
r(s)ds 



Gi- 



lt is assumed that such martingale exists, that is the set of given options satisfies the 

no-arbitrage condition. 

Remarks: 

(a) As before, we can assume r(s) = 0, since we can take hj — kije~ r ( s ) ds instead 



(x u -kj) + 



Cij. Also, we can take 



of kij , and rewrite all the equalities as E 
the function 

~g(X t .) := e-/o* ^ ds g ( e /o* ^X t .) 

instead of g(X t *). Notice that if g(X t *) is a continuous piecewise linear function, 
so is g(X u ). 

(b) Thus, again, only the order of t\, ti, ■ ■ ■ , t n is important, but not the specific values 
that the sequence takes. Thus we can assume {ti,t 2 , ■ ■ ■ , t n } = {1,2, ... ,n} for 
convenience. 

(c) Recall that we must have E[X t ] = X , where X is the current price of the stock. 
In light of the above remarks, Problem (b) can be reformulated as follows: 
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Problem (b) 



Given a continuous piecewise linear function g : R + — > M + and an ordered set 

S := {(fcy, C tj ) | (fc tJ -, C ti ) et^€{l,2,..., n}, je{l,2,..., U(t)}} (10) 
find the domain of values of 



E\g(X t .)] 



(11) 



for some t* € {1,2,..., n}, under the condition, that {Xl, . . . , X t * , . . . , X„} is a nonneg- 
ative martingale such that 



E 



(X t - k tJ ) + ] = C ti for all (k tj , C tj ) € S, 



(12) 



and 



£[X t ]=X for alH= l,...,n. 



3.1 Treatment of future conditions 

Let 

J 7 = {(fctj, C^) |i* < t < n, (k t j, Ctj) is a vertex of the convex hull of <S5° t , } 

represent future points which are on the border of the convex hull of 5>t* ■ And let 
(kfutj, Cfutj), j — 1, 2, ... , U(fut), be the elements of T . Then we have the following 
lemma: 

Lemma 1. If we reduce conditions (fT2"|) in Problem (b) to the set of conditions 



E 



E 



(Xt-ktjV =C tj , t<t*, j= 1,2..., U(t), (13) 



(X t , - k fut ,j) + <C futJ , j = l,2,...,U( fut ), (14) 



the solution to Problem (b) will not change. 

Proof. Suppose that we want to maximize E[g(X t *)] and g\ is the optimal solution 
under conditions (fl2)) and <?| i s the optimal solution under conditions (fl3)) and (fl4j) . 
Then g* < gj- L et us show that g{> g^- 



Let LT* be the distribution that maximizes J^°° g(x)dH(x) over distributions II of X, 



under conditions (fT3"|) fc tH}. Then there exists a martingale X\,... ,X t * which satisfies 
the part of conditions (fT2f for t < t* and with the distribution LT* of X t * , such that 
^n* (kfutj) < C(kf u tj) for all j — 1, . . . , U(f ut ). Since \&n* has the properties stated in 
Proposition[TJ none of the points (k, ^n* (k)), k > 0, is in the interior of the convex hull of 
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SS? t , and *&(k) < Xq — k. Then, since also the set S>t* by itself satisfies the no-arbitrage 
condition, by Theorem [4] there exists a martingale X? , . . . , X tn which satisfies the set 
{S>f }U{(fc, "fri* (k))} k>0 - Notice that the distribution of X t * is! ^n*- As in the proof of 
Proposition!!! it can be shown that the martingale X\, X2, ■ ■ ■ , X t * can be extended to the 
martingale X\, . . . , X-t* , . . . , X n by specifying transition probabilities C(Xt*+i\Xp) = 
C(Xf+i\Xt*), ■ ■ ■ ,Z{X n \X n -x) = C(X n \X n -i). This extended martingale satisfies all 
the conditions (12]). Thus g\ > g\. □ 

3.2 Restriction to discrete distributions 

From now on we make an additional assumption that 

P(Xi<L,i=l,...,t*) = l. 

for some L > 0. Let 

K s = {0} U {L} U {k tj I (fc y , C tj ) £5,(<f}U {k t] I {k tj , C tj ) e T] 

and 

K g = {ke [0,L] \g'(k-)^g'(k+)} 
be the set of points, where g changes its derivative. Define: 

K, = K-s U IC g . 

Theorem 5. Let X\, . . . ,Xf be a nonnegative Markov martingale with support in 
[0, L\. Then there exists a martingale X\, . . . , X t * with support in /C, such that 



E 



{x t -k t3 y\ =e 

E\(X t *-k futd ) 



4 = 1,2,. 



= E 



X»-k 



fut,j 



.,t*; j = l,2,...,f/(4) 
j = l,2,...,C/( /nt ) 



and 



E[g{X t ,)] =E[g(X t ,) . 
We will use the following two lemmas to prove Theorem 

Lemma 2. Let X\, X2 be a martingale with values in M. d . Then there exists a continuous 

time martingale I X t > with continuous paths, such that X\ = X\ and X2 = Xo. 
y > te[i,2] 

Proof. By the martingale version of Skorohod embedding theorem [131 p. 229], there 
exists a Brownian motion B and optional time t > 0, such that X2 — X\ = B T a.s. 
Let M t — X\ + BtAr- Then by optional stopping theorem M t is a martingale. This 
martingale is bounded and, thus, uniformly integrable. Consequently, there exists 
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such that {M t } tg j +oc j is a martingale and M t converges to almost surely by Doob's 
theorem O p.285]. 

Let h : [0, +00] — * [1, 2] be a nondecreasing continuous function and let M t = A^ft(t)> 
t G [1,2]. Then {Mi,M 2 } has the same joint distribution as {Xi,X 2 } and {Mt} te ^ 2 i 
is a continuous time martingale with continuous paths. □ 

Lemma 3. Let {-Xt} te n 2 ] be a continuous time martingale with values in M. d . Let 
ft : M. d — > R, be a continuous piecewise linear function. Let R d = Uj-Dj, and ft, is affine 
on each A, i = 1,2, ... ,n. Then for each i 



E 



{*t} t6[1)2] e A = -B ft(* 2 )| {^} te[1>2] e D. 



Proof. The result immediately follows from conditional Jensen's inequality applied to ft 
and —ft restricted to Di j9] p. 277]. □ 

Now let us return to the proof of Theorem 

Proof. We start at time t* and proceed recursively in the following way. Let B t , t G 
[0, +00), be a Brownian motion starting at 0. Let X t = X t * + B t _ t * for all t > t*. 
Then X\, X2, ■ ■ ■ , X t * , Xt is a martingale. If ti is the first hitting time of X t on /C, then 
by the optional stopping theorem X\,X2, ■ ■ ■ , Xf , .XW«'+i)ATi , X( 4 * +2 ) Ari , ■ ■ • is also a 
martingale. Since this martingale is bounded, it is uniformly integrable, and thus 



p. 283, p.285] there exists Xoo : such that E 



Xno On 



= X„ for all n = t*,t* + 1, . . ., and 



X„ converges to a.s. Let X t * = X, x . Then Xi,X 2 ,.. . ,Xt*-i,Xp is a martingale 
and P (Xr G /c) = 1. 

Now let us make the discrete time martingale Xi, X2, . . . , X t »_i, Xf be continuous on 
[t* — 1, t*] and have continuous paths. That is, by Lemma[2]we can define the martingale 
such that Xf-t — X t *-i, X t * — X t * and it is continuous. Let r 2 be the 



{*}, 



te[i*— 1,4*] 

stopping time which stops the martingale when it first hits K, after t* — 1. Then r 2 is 
bounded and Xt*_i, X T2 is a martingale by the optional stopping theorem. By Lemma[3] 
Xf-i (= X t *-i^ and X T2 takes the same values on functional but P (x Tt ,_ ± G K) = 1. 
Since Xi,X 2 , . . . , Xt*_i, X T2 , Xf is a martingale, so is Xi,X 2 , ■ ■ ■ ,Xf- 2l X T2 ,Xf. We 
let Xf — i = X T2 and proceed in this way until time 1. □ 

We next introduce Problem (b') that is equivalent to Problem (b). 
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Problem (b') 



Find the domain of values of 



under the condition, that {X\, . . . , X t *} is a nonnegative martingale with support in 
K := ICs U IC g such that 



E 



(X t -k tj ) + =C tj , t <t*, j = 1,2..., U(t), 



(15) 



E (Xf - kf ut .j) + <Cf ut< j, j = l,2,...,U(f ut ). 



From Lemma [T] and Theorem [5] we obtain 

Corollary 2. If the price of the stock is bounded by L, as in (|15p. then Problem (b) 
and Problem (b') have the same solutions. 

Proposition 5. If there is no bound on the price of the stock, then as L — > oo, the 
optimal solution of Problem (b') converges monotonically to the optimal solution of 
Problem (b). 

Proof. Notice, that since g is a continuous piecewise-linear function, E [g(Xf)] is always 
finite once E[X t *] is finite. Thus, the best possible upper bound on E [g(X t *)], C*, 
for Problem (b) is finite . If C* L is the best possible upper bound on E[g(Xf)] for 
Problem (b'), then C* L < C* and C* L approaches monotonically C* as L — > oo according 
to Corollary [21 Thus C* L — * C* as L — > oo. The same is true for the lower bound. □ 

3.3 Algorithm for the one dimensional case 

We have shown that in order to solve problem (b), it is enough to consider Markov 
martingales with support in K,. The distribution of a discrete time Markov process is 
fully determined by joint distributions of each two consecutive states. 

Let K. = {ki, &2, • • • , k n , . . . , k^} and K denote the set of all the indexes 1,2, ... ,N. 
Then the distribution of the Markov process with support in K. is determined by the 
collection of 

ftim.na =P(X t = k ni ,X t+ i = k n2 ), t = l,...,f, ni,n 2 eN. 

Since most of the constraints that we consider are constraints on marginal distribu- 
tions of the process, for convenience we also introduce 

Pf,n = P{X t = k n ) t = l,...,t\ neN. 

For P to be a measure, it must satisfy the countable additivity condition: 
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E p t;n, ni = Pt+i;m for all t < t* , n x G X, 
new 

and nonnegativity condition: 

ft ; n>0, P t:ni ,n 2 >0, for all t <t* , n, tii, n 2 G X. 
In addition, for the probability measure we must have 

E P Un = I- 

neN 

The martingale condition can be written as 

E P t;ni,n fc n = ^t;ni*m for all * < <*j n l G X. 
new 

For the martingale to satisfy prices of options with maturities less than or equal to t* , 
we must have: 

E P ';n(fcn - h )+ = C tj for all t < t*, j = 1,2..., U(t) 

neW 

As was proved in Theorem [5l for prices of options of maturities greater than t* to be 
satisfied, it is enough to have 

E p t* Akn - kf u t,j) + < Cf u t :j for all j = 1,2 ... , U( fut ). 
new 

Thus the problem of determining best possible upper (lower) bounds on the price of 
an option with maturity t* and a payoff function g, consists of finding the collection of 
Pt-n and Pt ; ni,n 2 i such that 



max(min) ^ Pt* ;n g(k n ) 
new 

s.t. = X > 

neW 

E P *;m,n = Ptjni for all t < t*, m G X, 

E P *;n,ni = p *+i;ni for all i < £*, n x G X, 

new 

E p t;n u nk n = Pp,mk ni for all i < **, tlx G X, 
new 

E - M + = C tj for all * < **, j = 1, 2 . . . , U(t), 

new 

E • P t*;n(fcn - kf u t,j) + < Cf ut j for all j = 1, 2 . . . ,U( fut ), 



new 

-Pt;ni,n 2 > for all f < til, tl2 G X. 

This is a linear optimization problem with 0(N 2 ), where iV is the number of elements 
of K. 
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4 The Two Dimensional Case 



In this section, we address Problem (c), i.e., the problem of determining bounds on 

options on two different assets if prices of options on individual assets are known. 

Let \ S (t) \ be a stochastic process describing the price process of the first 
I J t>o 

asset, and |5"(i)j that of the second asset. If (fc^-, C^, ij), i = 1, ...,n, j 

1.2 U(t, h), represent the set of options on asset S h , with ti being the maturity of 

an option, its strike price, and its price, then we must have 



C 1 



"(s)ds 



Ec 



( 



X\\ - khe 



-(s)ds 



\ 



, (16) 



where {X^] t>Q , is a martingale under Q. As in the single dimensional case, we can 
assume without loss of generality that the risk-free rate r(s) can be assumed 0. Thus, 
we can change times t% to their indexes. 

LetT := {(h,t,j)\h G {r,n }, i G {1,2,.. .,n},j = 1,2,. ..U(t,h)}. We define a fam- 
ily of measurable functions : — * M+ as 

/^.(X t ) = (X t fc -^.) + . (h,t,j)ET. (17) 
Then the problem can be formulated as follows: 

Problem (c) 

Given the set (fcy,C|l-), (h,t,j) £ T, and a continuous piecewise linear function g : 
R 2 — > R, find the maximal and minimal possible values of E[g(X.f>)], for a given t* G 
{1, 2, . . . , n}, under the condition that Xi, X2, . . . , X„ is a two dimensional nonnegative 
martingale, such that 



E[f t )(X t )} =C%, {h,t,j)eT 



(18) 



Definition 3. We call function g : M. d — * R piecewise-linear, if there exists a partition of 
R d by a finite number of hyperplanes (lines if d = 2) into nonoverlapping subsets (regions 
if d = 2), such that in the interior of each subset g is an affine function. 

Notice, that by Proposition [4] we might restrict our consideration to Markov martin- 
gales only. 

4.1 Treatment of future conditions 

Notice that the function ffy depends only on the /i-coordinate of X t . With this property 
of /—constraints, we can show that future conditions have a "simple effect", in the sense 
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that conditions (fl8|) for t > t* can be reduced to separate conditions on the first and 
second coordinates of X t » . Let 



S£ t . = { C*) |t* < t < ra,t e {1,2,..., C/(t, ft)}} , ft , 
5^S = 5| t . U {(0, +00), (+00, 0)}, ft =, ,u , 

and 

J 17 ' 1 = { (k' t ),C^) \t* <t< n, (k^Cl)) is a vertex of the convex hull of } 

represent future points which are on the border of the convex hull of 6>>^? . Let 

(kfut,j>Cf ut j) , j — 1, 2, ... , U(f u t, ft) 

be enumerated elements of T h . 

In formulating all further results we will assume that \k t pC t ^j, 

(i,t,j) <G T and ffcy,CyJ, (u,t,j) G T separately do satisfy the no-arbitrage con- 
dition, that is that there exist one-dimensional martingales X t and X t such that 



(X' t l -k' t )) + =C' t ), h=!,U 



Notice that this condition ensures that the two dimensional martingale which satisfies 
conditions (1181) also exists. 



Theorem 6. Let Xi, X2, . . . , Xi* be a two-dimensional martingale satisfying 
E [/tj-(X t )] = C£j, for all {(h,t,j) G T | i < t*}. If also for each ft G {1,11 } 



E 



{X?,-k h fut ^]<C h fuK 



is satisfied for all (kf ut a , Ch ut j ) G JF' 1 , then there exist two-dimensional random variables 
X t » + i, . . . , X„ such that 

(a) Xi,X 2 , . . . ,X t .,X t . + i, . . . ,X„ is a martingale; 

(b) all the conditions E [/A(X t )] = Ch are satisfied for all (h,t,j) G T. 

Proof. As in the proof of Lemma [TJ we can show that for each ft G {1 ,n } there exists a 
one-dimensional martingale X/l , Xji +1 , . . . , X^ 1 , such that it satisfies the set (fc^-, C t ^) , 
i > **, J = 1, . . . , U(t, ft) and * ^ h = W X h . 

We will extend the martingale Xi, X2, . . . , X t * to the martingale Xi, X2, . . . , X t * , 
Xf+i, . . . , X„, such that X t * , ~K t *+i, ■ ■ ■ , X„ is a Markov process. To define a joint 
distribution of X 4 and X i+ i for t > t* , allow each coordinate to evolve independently 
of each other. That is, for each ft =1,11 define C (^X^\X^ +1 ) , as C (^X^\X^ + -^j for t — 
t*,t* + l,...,n. □ 
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Theorem [6] allows us to reduce "future conditions" to conditions at time t* . This is 
a significant simplification of the problem, since difficulties arise from the requirement 
that the martingale satisfies E [//j(X t )] = Ck, (h,t,j) e T. 

Corollary 3. Our main problem is equivalent to the following: find the maximum (min- 
imum) of E [g (X(«)] under the condition that Xi, X2, . . . , X t * is a Markov Martingale 
such that 



(a) E 

(b) E 



= ty, for all t<t*,j< U{t, h) and 
" & M j) + ] < C)ut,i for all (k h fut>j ,Cl t>j ) G T\ 



Definition 4. Let f) ut% . : K^. — > R + , /) utJ (X) = - A^J , be a family of 
functions, where each function corresponds to some k^ ut j, such that (j$f u t j>Cf u t j) e 

Constraints from the past introduce much more difficulty to the problem than fu- 
ture constraints, ft is not possible to reduce past martingale conditions to separate 
conditions on X t , and X t , . In other words, for a given distribution of X t « , and two 
sequences of marginal laws {n f } 4=1 .„ and {Qt} t= i t « of one dimensional processes, 
such that Il t « = C fe) and Q t , =£ (Xt'^j , the existence of one-dimensional martin- 
gales with the above sequences of marginal laws, is not sufficient for the existence of a 
two-dimensional martingale Xi, X2, . . . , Xt* such that C (^X t ^ = lit and C (x t j = Qt 
fori = l,...,i*. 

Counterexample. Suppose that X2 has an atomic distributions with two atoms a and 
b, as shown on Figure [SJ where each atom has weight of 1/2. Suppose that the marginal 
law of X 1 is P (Xi = x^j = 1, where x c — ^(x a + Xb) and the marginal law of X 1 

is P (^X'i = ya^j = P (X-± = ybj = \- Obviously, we can define a joint distribution 
of X 1 , X 2 , so that it is a one-dimensional martingale. The same is true for X 1 , X 2 ■ 
However, it is not possible to define a joint distribution of Xi, X2, so that it is a two- 
dimensional martingale. 

Thus, the condition that Xi, X2, . . . , Xt* is a martingale, is truly a "joint" condition 
on {X' t } and{x t IJ }. 

4.2 Restriction to discrete distributions 

In what follows, function g is defined as g : — *■ R + , (x,y) — > \ax + fly — k g \ + . 
However, an analogous solution would apply if g : — > R is any continuous piecewise- 
linear function. 

As in the one dimensional case we add an additional constraint that stock prices are 
bounded, that is, there exists L > 0, such that for h = I, II 

P(X' 1 \X' 2 \... 1 X':<L) =1. 
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Vb 

Va 

•&a %b 
Figure 5: Counterexample. 

Definition 5. For each t = 1, 2, . . . , t* — 1 we define Graph Gt (Figure [6]), which is 
contained in [0, L] x [0, L] and has vertices on intersections of the lines: 
(a) x = k\ v j < U(t,i ), (b) y = j < U(t,u ), 

(c) segments of Q: from (0,0) to (0,L), (0, L) to (L,L), (L, L) to (L,0), (0,0) to (L,Q). 

For t* we define Graph Gt* (Figure [7]) which has vertices on intersections of the lines 

(a) x = k'fj, j < U(t,i ), (a') x = k' futd , j < U( fuU i ), 

(b) y = kll p j < U(t,u ), (b') y = k'/ utJ , j < U( fut ,u ), 

(c) segments of Q: from (0,0) to (0,L), (0, L) to (L,L), (L, L) to (L,Q), (0,0) to (L,Q), 

(d) ax + (3y — k g . 

Each edge of Gt belongs to one of the corresponding lines (or segments). The inter- 
section between two edges is either empty or consists of one vertex point. 

k"i k"k 




fc 7 ^ ak'-ffk"'---.. k' 



Figure 6: The form of Graph G t for t = Figure 7: The form of Graph G t * ■ 

l,2,...,t*-l. 

Let 

• G\ designate all the vertices of Graph G t , 
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• Gf denote the set of points which lie on edges of Gt and 

• Gf- denote regions of K 2 over which edges of Graph Gt splits [0,L] x [0, L\. 

Let us enumerate regions of Gf . We will thus write r for the region of Gf numbered 
r. Let 3? t be the set of all integers such that each of them represents some region of Gf 
and V be the set of all integers such that each of them represents some vertex of GjC . Let 
5Heg t : [0, L] x [0, L] — ► Ji t , be the function which attributes to each point of [0, L\ x [0, L] 
the index of the region that contains it, and let Ver : [0, L] x [0, L] — > V, be the function 
which attributes to each vertex of GjC its index v. Since edges are contained by several 
regions, we have to specify which region corresponds to each edge. 

Definition 13. We will say that a random process Xi, X2, . . . , X 4 * has property Pq if 
and only if all of the following hold 

(a) P(X t . eG£) =1; 

(b) Xi has an atomic distribution with at most one atom per region of Gf ; 

(c) If n E Dli, r 2 S ft 2 • • ■ , n e % with t < t*, then 

£ (X f+1 1 pReg^Xx), <Ke flt (X 2 ), . . . , 9tefc(X,)] = [n, r 2 , . . . , r,]) 
has an atomic measure with at most one atom per region of Gf . 

Theorem 7. Let Xi, X 2 , . . . , X t » be a nonnegative Markov Martingale bounded by L. 
Then there exists a nonnegative Markov Martingale {Xi, X2, . . . , Xt* }, bounded by L, 
such that it has property Pq and satisfies 

e [/£.(x t )] - e [/^ (x t )] ,*<**, (M,i) e 

^[/^■(X*,)] =S[/ / \ tJ (X t .)], i = l,2,...,?7( /ut) /i), 

S[ S (X t .)]=S[ 9 (x t .)]. 

Proof. For i = 1, . . . , £*, let (Tj designate the smallest a- algebra for which Xi, X 2 , . . . , X,; 
are all measurable. 

The proof goes by induction. Let a\ be the smallest er-algebra for which the ran- 
dom variable 5Heg 1 (Xi) is measurable. Then o\ C <j\. Define Xi = ^[Xila-i]. Then 
{Xi, Xi; CTx> cr il is a martingale. Xi has an atomic distribution with at most one atom 
Zi- r = E [Xi \ (Xi) = r] per each region r of Gf and takes the same values on 
functionals /•&, j = 1, 2, . . . , 17(1, ft.), ft G }, as Xi by Lemma[3] 

Now, Xi,X2, . . . ,Xf is also a martingale. Let tr 2 be the smallest er-algebra which 
includes o\ and for which *Heg 2 (X 2 ) is measurable. Define X 2 = i?[X 2 |cr 2 ]. Then 
{Xi, X 2 , X 2 ; cti, cr 2 , cr 2 } is a martingale. Indeed, E [X 2 |cti] = E [E [X2|cr 2 ] \a\] — 
E [X 2 |cti] = Xi, since a\ C a 2 , and E [X 2 |ct 2 ] = X 2 by definition. C (X 2 ffitQi (Xi) = r ) 
has an atomic measure with at most one atom Z 2;r%f = E [X 2 I^Kefli (Xi) = r, 9^eg 2 (X 2 ) = f ] 
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per each region f of G^. Also by Lemma [3] X2 takes the same values on the functionals 
j = l,2,...,C/(l,ft), ft G {/,«}, as X 2 . 

Continuing in the same way, we can form a Markov Martingale Xi, . . . , X t *-i, X t » 
such that £ (X t | [9 f le0 1 (Xi), (X 2 ), . . . ,9tefl t _ 1 (X 4 _i)] = [ri,r 2 , . . . ,r t -{\) has an atomic 
measure with at most one atom per region of Jit for all i = 1,2, ... ,t*. 

Now, let Bj_t« be a two-dimensional Brownian motion starting at 0. We define 
for t > t* X t = X t * + B t _ t . and let r be the first hitting time of X t on Gf, . Then 
X t », X t * Ar , X( t * +1 ) Ar , . . ., is a martingale by the optional stopping theorem. This mar- 
tingale is bounded, thus there exists Xqo, such that X( f » +1 ) AT , X( t * +2 ) Ar , . . . , Xoo is a 
martingale and X( f » +1 j AT , X( f » +2 ) AT , . . . , converges to X^ a.s. by Doob's theorem. Let 
X(» = Xqo. 

By Lemma H P (x t * G Gf,^j = 1 and X t . takes the same values on functionals 
E[fp,j{.)], E[fj ut -(.)] and E[g(.)] as X t *. In the same way, but by adding an one 
dimensional Brownian motion to X t * over the edge to which X t * belongs, we can reduce 
the support of the martingale at time t* to Gj* . □ 

Thus, to solve optimization problem (c), it is enough to consider only Markov martin- 
gales Xi, X 2 , . . . , X t . satisfying the property Pq- Consequently, it is a finite dimensional 
optimization problem (distributions with the Pq property can be described with finitely 
many parameters.) General measures belong to an infinite-dimensional space which is 
even uncountable. So the above result is a significant simplification of the problem. 
However we will see that the number of parameters grows exponentially with t* . So 
solution methods described below will work only if t* is not very large. 

4.3 Algorithm for the two dimensional case 

By Theorem [7] in order to solve Problem (c) it is enough to consider martingales with 
property Pq. That means that the distribution of the initial state (i = 1) is atomic 
with at most one atom per each region of Gf , the distribution of the final state (t = t*) 
is atomic with atoms at the vertices of G t » , and the conditional distribution of each 
intermediate state (1 < t < t*) is atomic with at most one atom per each region of Gf ', 
conditioned that the sequence of all prior visited regions is known. 

The distribution of the described process is fully determined by the probabilities 

^n,..,^-!,^ = P ([^ e fli(Xi)i • • . ,5tefl t ._ 1 (X t ._i), Ver(Xt.)] = [n, . . .,rf-i,v]) i 1 ®) 

of visiting all possible sequences of regions (vertices at the final state), and by the unique 
places of visits of each particular region, conditioned on the prior sequence of regions. 
Let Zt- ri ,r 2 n G denote the place of the unique atom of the conditional distribution 

£(X t I [«efl 1 (X 1 ) > ... > 5tefl t (X t )] = [ri,r 2 ,...,r t ]) > 
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where r\ £ %x, r 2 <E 9?2, ■ ■ ■ ,r t e R t , t < t* . Since all regions of % for t = 1, 2, . . . , t* — 1 
are rectangular, Zt. ri r2] ... irt can be represented as 



Z{;n,r2,...,r i — u r t + ' w r t z t;r 1 ,r 2 ,...,r t v r t z t;r 1 ,r 2 ,...,r t i 



(20) 



where u r( , wy t , v rt £ are vectors uniquely defined by the rectangular rt (see Figure 



D and 



e [0, l], /i =j,j 



A" 



k' 

Figure 8: Parametrization of a rectangular. 

Notice that / tJ , / tJ are affine when restricted to one of the regions of Gf, and 

/& (Z t . ri) ..., rt _ lir ) = a t ^(r)4 ri!r2) ... irt _ l!r + 6* (r), r G ft* (21) 

for some coefficients aj^ (r), &y(r), which are easy to determine. 
For convenience we also introduce 

l'>, . ,. - P ([^(Xi), . . . ,9te flt (X t )] = [n, . . . ,r t ]) 

for each t < t* , which are probabilities of visiting a particular sequence of regions prior 
to time t. 

In order to satisfy the conditions of a probability measure, in addition to nonnega- 
tivity, we must have 



(22) 



Pt— l;ri,r2,...,i"t_i ^ ] Pt;r 1} r 2 ,...,r t - 
relit 



ur , for all n € 31i,...,r t _i €# t -i, * = 2, . . . , t* - 1, 



Pt»-l;r 1 ,r- 2 ,...,r t »_ 1 = ^ Pt i r i. r !-. r t*-i,». for all 1"i € . . . , E Rf-l, v € V. (24) 



To satisfy the martingale condition we must have 

Zt_X;nr2,...,rt_i-ft— lin,n,..,,rt-i = X/reX Zt;rir 2 ,...,r t _i,r-ft; 

for all n e 3?i, . . . ,r t _i € 3? f -i, t = 2, 3, . . . , t* - 1, 



(25) 



2G 







1 



z 



Figure 9: Linear inequalities on 2t;ri,r 2 ,...,r t _i,r t and P< 



t;r\,T2,.~,rt ■ 



Zt*-l ; rir2,...,r t »_ 1 -ft— l;ri,r 2 ,...,r t *_ 1 — SueV Zi*;r 1 r2,...,r t *_i 1 i)-Pt*;ri,r 3 ,...,n-i>'i> (26) 

for all T\ eS b .. . ,r t ._i € w G V. 



Notice, that Zt ;rir . 2! ... jrt _ 1)7 . is defined in terms of z\. r r2 r r and z" r r r r in (|20)l . 

Finally, the measure must satisfy the prices of options. We use the representation (|2"Tj) 
of payoff functions to write down the constraints from options prices. For satisfying past 
constraints we must have 

SriG3?i,...,r t eK t ( a tj ( r t ) z t\n ,r 2 ,. . . ,rt- 1 ,rt ( r <)) Pt;r l ,r 2 ,...,r t = j (27) 

for {(h,t,j) ET\t<t*.} 

The present constraint are 

SnG3Ji,...,r t ._ 1 63l t *_i,«eV ffj( V )Pf;r x ,— ,r t * _ x ,v = Cfji (28) 

for {(h,t,j) &T\t = t*.} 
Finally, constraints from options of maturities larger than t* are 

SrieKi,...,r t »_ 1 e3i t «_i,«eV ffutj(. v )Pt*;r 1 ,r 2 ,...,r t ,_ 1 ,v < @fut,j ft®) 

for {(/i,i,j)GT|t>r.} 
Thus problem (c) can be written as : 

Find the collection of Pt-. ri ,r 2 ,...,r t > and zb u €[0,1], ne r t €%, t = 

1, 2, . . . ,t* — 1, /i =i,n, such that it satisfies conditions (f22 |) - ([29|) and 



rie» 1 ,...,r t ._ 1 e3J t »_ 1 ,tieV 



is maximized (minimized.) 
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The above system is not linear. By making the change of variables 



z t;r 1 ,r 2 ,... : rt-i,r t ~ z t;r 1 ,r 2 ,...,r t -i,rf^t; 



r 1 ,r 2 ,...,r t 



(30) 



for t = 1, 2, ...,£* — 1, we make it linear. The domain for {zpri r2,...,r t t Pt;ri,n ...rt) * s 
represented on Figure [5] Thus all the constraints as well as the objective function are 
linear. 

4.4 Computational Example 

The data on VerizonCM and Cisco call options (Table [J) is taken from The Wall Street 
Journal of October 22 2002 to derive best possible bounds on the price of an option with 
a payoff function (V + C — k) + , where V is the price of the VerizonCM stock at the 
maturity, C is the price of Cisco stock and k is the strike price. The prices of VerizonCM 
and Cisco stocks on that day were $37.75 and $11.22 correspondingly. 

Table [2] represents best possible bounds on the price of the option with maturity 52, 
if all given options are taken into account. Table [3] gives bounds on the price of the same 
option, if only options of maturity 52 are taken into account. As expected, these bounds 
are looser. 

As can be noticed, bounds for strike prices smaller than 20 are extremely tight. That 
makes sense, since the probability that the sum of prices of the two given stocks will go 
below 20, is almost zero. Thus, the option will be exercised almost surely and buying 
the option is almost equivalent to buying the two stocks. 

4.5 Approximation approach 

In this section, we develop an approximation approach which allows to overcome the 
exponential growth (in t*) of the number of variables. 

Recall that we consider only martingales X t such that P(X f < L) = 1 for some L > 0. 
For now we will allow g : K.+ — > K., to be any function continuous on [0, L] x [0, L] and 
develop a general approximation algorithm. Then we will show how the efficiency of the 
algorithm can be improved, if g is a continuous piecewise linear function. Choose e > 
and designate by L e the set of all vertices of the e-square lattice in [0, L] x [0, L) Cl 2 . 

Theorem 8. Let Xi, X2, . . . , X f » be a two-dimensional Markov Martingale with values 
in [0, L] x [0, L] . Then there exists a two-dimensional Markov Martingale Xi , X2, . . . , Xf 
such that for alii = 1, ...,**, P (x t 6l' E )=l and 




h 



=1,11 



Here L\ 



is an e-square lattice in [0, L + (t — l)e] x [0, L + (t — l)e]. 
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VERIZON 


CISCO 


maturity 


strike price 


price 


maturity 


strike price 


price 


17 


30.0 


7.80 


17 


10.0 


1.65 


17 


35.0 


3.70 


17 


12.5 


0.35 


17 


37.5 


2.05 


52 


12.5 


0.65 


17 


40.0 


0.95 


80 


7.50 


3.75 


17 


42.5 


0.35 


80 


10.0 


2.15 


52 


40.0 


1.80 


80 


12.5 


0.80 


80 


30.0 


8.60 


80 


15.0 


0.30 


80 


35.0 


5.10 


171 


15.0 


0.65 


80 


40.0 


2.20 








171 


42.5 


2.00 









Table 1: Data. 



strike price 
lower bound 
upper bound 


10 20 30 40 45 50 55 60 70 
48.97 38.98 28.99 19.09 9.02 4.03 0.89 0.57 0.04 0.01 
48.97 39.02 29.31 19.47 10.56 7.17 3.91 2.08 1.93 1.62 


Table 2: Best possible lower and upper bounds on the price of an option on the linear combi- 
nation of two stocks. Maturity is 52. Options of all maturities are taken into account. 


strike price 
lower bound 
upper bound 


10 20 30 40 45 50 55 60 70 
48.97 38.98 28.97 18.97 8.97 3.97 0.25 0.00 0.00 0.00 
48.97 39.97 30.99 21.99 13.01 8.79 4.56 2.37 2.23 1.92 


Table 3: Best possible lower and upper bounds on the price of an option on the linear combi- 
nation of two stocks. Maturity is 52. Only options of this maturity are taken into account. 
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Proof. Let B\ C B2 C • • • C Bt* be the family of cr-algebras corresponding to the 
martingale Xi, X2, . . . , X t » . Let A t := X t — X t „i for t = 2, 3, . . . , t*. Then 



X t = Xi + A 2 + A 3 



andf;[A t |B t _i] =0. 





e 

L ^ 


SiW 




Six) 










Y 










Six) 




S 4 (x) 














































-L 



Figure 10: Lattice L e 



Let L £ designate the set of all vertices of the e— square lattice in [— L, L] x [— £, L] G K 2 
(Figure [TOj) And let : [— L, L] x [-L, L] L e , j = 1,2,3,4, be the family of 
functions, such that each function puts in correspondence to x G [—L,L] x [—L,L] one 
of the closest to it four vertices of L e . 

Notice that A t G [—L,L] x [— i,L] for all t = 2, ...,£*. Let us define a random 
variable At such that conditioned on A t , it can take only four values: {5 J '(At) — A 4 , 



j = 1, 2, 3, 4} and such that -E 



A t 



B, 



= 0. Moreover, we add a requirement, that con- 
ditioned on Bt, At is independent of any other random variables in consideration, includ- 
ing Ai, A 2 , . . . , At-i. While A 2 , A 3 , . . . , A t * can be defined based on A 2 , A 3 , . . . , A t * 
as described above, define Ai in a similar way, but take Ai = Xi. 



Then we have A t + A t G L e and P 



(|A?| < e) 



— 1, h — t = 1,2. 



,t*. Let 



a ( A l5 A2, • • • , A t ) denote the smallest er-algebra for which Ai,A 2 ,...,A t are all 

measurable and 



Bt:=a B 4 U<7 Ai,A a ,...,A t ),t = l,2 



Xi =Xi + Ai, 

X 2 = X 2 + Ai + A 2 , 



X t . = X*. + Ax 



A t , 
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Xi € L e is obvious. Now 
X 2 = X 2 + Ai + A 2 = Xi + (X a - XO + A x + A 2 = ( X x + A x 



A 2 e L t 



But since also P i 



^2 " ^2 f 



< 2e 



1, /i 



and X2 £ £ £ C L e , it must be that 



X 2 € E\. Continuing recursively the theorem follows. 



□ 



If Xi, X 2 , . . . , X n is a two-dimensional mar tingale, that satisfies E [/^-(X t )] = C[) 
for (t,j, h) € T , and Xi, X 2 , . . . , X n is a two-dimensional martingale, such that 



P 



< te J = 1, h =1 ,n , t = 1, . . . , n, then we must have 



E 



m x 



< et, h =1,11 , t = 1, ... ,n 



and 



^[ff(X t .)] 



< My/let* 



<?(X t . 

where M is the smallest number, such that 

|ff(x 2 ) - ff(xj.)| < M||x x -X2II2 

for all xi, X2 > 0, 

Notice that if y) = |cke + /3y — fc| + , then \[2M = a + [3. 



(31) 



(32) 



Theorem 9. The following problem has its optimal solution converging to the optimal 
solution of problem (fP8|) as e — > 0: 

over all Markov martingales X 1; X 2 , . . . , X t * , such 



Minimize (maximize) E g ^X t *^ 
that P (X t € Llj = 1 and 



/A ( X t )1 - < et, {(h, t,j) eT\t< t*} , 



(x t . ) ] - c^ tJ < et* , { , i) e r 1 {k% tjj , c h futtj )er h ) 



(33) 



(34) 



Proof. Let us prove the theorem for the case when we want to maximize £%(X t »)]. For 
the minimization problem the proof is similar. 

Let g* be the optimal solution to the original problem (c) and g s the supremum as 
e — > of optimal solutions to problems with martingales' state spaces restricted to L\ 
and subject to the constraints (|33|) and (f34|) . From Theorem [8] and ([32]) it immediately 
follows that g s > g* . Suppose that there exists 5 > 0, such that g s > g* + S. Then there 
exists a sequence of e's, {{e n } | ei > £2 > ■ • • £n > ■ • • ; hnin-i-oo e n — > 0}, such that the 



limit of optimal solutions to problems corresponding to ei, e%, 



is equal to g s 
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Let {Xj™} n>1 be the sequence of martingales corresponding to the sequence {e„}, such 
that {Xf"} n>1 defines an optimal solution to the e„-problem. 
By definition of {X^™}, it takes values in L\ n and 

|£[/y(X t e «)]-Cy|<e n t, *=l,2,...,t*. 

Let //" designate the law of {Xj"} t=1 t , . Then /x £l ,/j, £2 , . . . , /i £ „, . . . have support in 
[0, L + t*ei} 2 , which is a compact. It follows that the sequence of laws {A t e rl }„> 1 is 
uniformly tight. Consequently, there exists a subsequence At £ll(c) — > (J- for some law /i [9, 
p. 230]. But then we have an admissible solution to our original problem with an optimal 
value greater than g* + S, which is a contradiction. □ 

4.6 Algorithm for the approximation approach 

Let us enumerate all the nodes of L\ and denote by 3ST the set of all the indexes. Designate 
by Z n € L\ , the n-th node of L\ . Let 

P t , n = P (X t = Z n ) and 

-Pt;ni,rt 2 — P = Z ni ; Xt+l = Z„ 2 ^ . 

Notice that the distribution of Xi,X 2 , ■ ■ ■ ,X t * is uniquely determined by Pt :ni! n 2 j ^ G 
{1,2,...,**}, m,n 2 eN. 

Following the same set of arguments as in the description of the algorithm for the 
exact solution to the one dimensional case, we can state that the problem consists in 
finding the collection of P f;n and Pt ;nil n 2 , t £ {1,2,...,**}, n, ni,ri2 G N, which solves 
the linear optimization problem: 

max (min) 2J P t*;n ff(Z tt ) 
new 

s -t- XmeW ^lai = 1) 

XI P «;ni,n = Pi ittl for all t < **, m € X, 
new 

X P *;n, ni = Pt+i;m for all t < t* , m e N, 
neW 

X p t; ni ,nZ„ = P t;ni Z ni for all * < f, m G X, 
new 

X P t; „/£(Z„) - < ie, X P ; „/£(Z„) - > -tc 
new new 

for {(/i,*,i)eT| *<**}, 

^ P t*;n//ut,j(Zn) - C/ut,j < 

new 

for {(M,j) G r I t > <*}, 

Pt;n!,n 2 > for all* < t* , til , n 2 £ N. 

As before this is a linear optimization problem. 



32 



4.7 Approximation approach for the special case of a payoff func- 
tion 



In this section, we show how the efficiency of the approximation algorithm can be im- 
proved, if g is a continuous piecewise linear function. 

Let G E — Gf U G E U ■ ■ • G E and let graph G be the graph with points of edges 
represented by G E . Then, if we restrict the state space of martingales to G E , the 
solution to problem (c) will not change. More precisely: 

Proposition 6. Let Xi, X2, . . . , X f » be a nonnegative two dimensional Markov Martin- 
gale bounded by L. Then there exists a nonnegative Markov Martingale |Xi,X2, ■ ■ • ,Xf|, 
with the state space in G E , such that it satisfies 



£[/fc(X t )]=J5[/£(Xt) 

E [fU,3 ( X *)]=£[/jW(Xt 

E[g(X t ,)]=E 



,*<**, (h,t,j)eT, 
j = l,2,...,U( fut ,t*), 



(x,.) 



Proof. The proof is similar to the proof of Theorem [5] except that we deal with two 
dimensional random variables and, thus, a two dimensional Brownian motion should be 
introduced instead of one dimensional. □ 



Let 



A* := {(x,y) e L\ | 3 (x ,y ) £ G E : |a?-x | < te, \y - y \ < te} . 

Then Theorem [9] will hold if we take A' instead of L\. The number of variables in 
this case grows only linear versus quadratically in the case of a general function g. 



4.8 Generalization to a multi-dimensional case 

The exact solution, the approximation approach for an arbitrary payoff function g and 
the approximation approach for a continuous piecewise linear payoff function g, can all 
be extended from a two dimensional case to a multiple dimensional case. All theorems 
and definitions of this chapter can be reformulated for an n-dimcnsional case, and proofs 
will be identical up to a dimensionality. However the number of variables in the linear 
optimization problem will grow exponentially in the dimension. 

4.9 Generalization to the case of options with continuous piece- 
wise linear payoff functions 

Across the paper we considered the case when prices of simple call options are given and 
bounds on prices of options with continuous piecewise linear functions must be found. 
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From financial point of view this formulation of the problem can be motivated by the fact 
that simple call and put options are more liquid and prices on them are readily available 
in the market, compared to more complex exotic options. Notice, that conditions on 
European put options can be easily expressed as conditions on call options. 

However, the problem could be extended to a more general case when options with 
continuous piecewise linear functions are given and bounds on similar type of options 
are to be determined. If there were no options of maturities larger than t* , then all that 
has to be changed in the solution is definitions of the set K, in the one dimensional case 
and graphs Gt in the two dimensional case. In particular, in the one dimensional case, in 
section 3.2, K, must include all points of ffi. + where payoff functions of given and target 
options change their derivatives. Similarly, for the two dimensional case, in section 4.2, 
lines defined in (a) and (b) of Definition [5] by strike prices of call options would have 
to be defined by points of M + where payoff functions of options on each asset change 
their derivatives. Then theorems that refer to definitions of set K and graphs G t would 
remain valid. 

However, the simplification that was made for treatment of future conditions (sections 
3.1 and 4.1) is not possible if payoff functions of options of maturities larger than t* 
are piecewise linear, since these functions are not necessarily convex. As a result, the 
martingale has to be constructed not only up to time t*, but up to the time of the 
maximum options maturity, and future and present constraints have to be handled in 
the same way as past constraints. In the one dimensional case, since the state space of the 
martingale at each time is the same, there is no principal difference. All that one has to 
do in addition to chang es in the definition of set /C is to extend the set of unknowns Pfn 
and -Pt:rti,n 2 up to the maximum time t and change inequalities to equalities in the future 
constraints. The problem remains a linear optimization problem. In the two dimensional 
case, however, there would not necessarily exist a martingale that takes values at the 
vertices of graph G t * at time t* and satisfy all future constraints. As was shown in the 
proof of Theorem [7j there does exist a martingale that satisfies all given constraints and 
at time t* has an atomic conditional distribution with at most one atom per each region 
of Gp , conditioned on the sequence of all prior visited regions. However, since regions 
of Gfi are not necessarily rectangular, the simple form of parametrization (|20|1 and the 
change of variables (|30|) can not be applied. Therefore, we don't think it is possible to 
reduce present constraints and the objective function to the linear form in this case. 

It is easy to see that the approximation approach for the two dimensional case can be 
applied to the solution of the problem, when prices of options with continuous piecewise 
linear payoff functions are given. In this case, one only needs to extend the set of 
unknowns Pt ;n and Pt ;ni ,n 2 from time t* to the maximum maturity and use equalities 
instead of inequalities for future constraints. 
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5 Static arbitrage bounds on basket option prices 

In this section, we address Problem (c) defined as follows: 

Problem (d) 

Given C G R™, k e R™, u, e 1", i = l,...,m and k > 0, w E R™, find the upper 
and lower bounds on 

E(^X-k ) + , (35) 

with respect to distributions of an n— dimensional random variable X with finite expec- 
tation and support in R" under the condition 

E (ujf X - h) + =Ci, i = 1, ... ,m. (36) 



To solve the problem we make an additional assumption that the support of X is 
bounded, that is, there exists L > 0, such that P [X h < LA = 1 for any h = 1, 2, . . . , n. 
Let fi : R™ — > R+ be the family of measurable functions defined as 

/ i (X) = (wfX-fc i ) + , i = 0,l,...,m. 

Let UjDj = [0, L] n be the partition of [0, L] n , into subsets Dj of R^L such that in the 
interior of each subset Dj all functions fi, i = 0,1, ... ,m are affinc. 

Let G C [0, L] n be the graph in R™ formed by intersecting hyperplanes 

w^X = hi, i = 1, 2, . . . ,m, 

X h = 0, h = 1,2,..., n, (37) 
X h = L, h=l,2,...,n. 

Let us enumerate these hyperplanes and call them hi, ■ ■ ■ , ^2n+m correspondingly. 
We call a vector re [0, L] n a vertex of G if there are n independent hyperplanes of (|3~T|) 
which intersect at r. Let Gy designate the set of all vertices. Notice, that this set 
contains no more than ff„i m = (2K+m)! e j emen ^ s _ 

Corollary 4. For each random variable X with support in [0, L] n there exists a random 
variable X such that P (X e G v ) = 1 and E [fi(X)\ = E [/<(X)] , i = 0, 1, . . . , m. 

Thus, to find optimal bounds we solve a linear optimization problem with unknown 
variables being weights of the atomic distribution. The atoms are located in the vertices 
of the defined graph. 
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